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A new method is presented here for evaluating approximately the pulsation modes of relativistic 
stellar models. This approximation relies on the fact that gravitational radiation influences these 
modes only on timescales that are much longer than the basic hydrodynamic timescale of the sys- 
tem. This makes it possible to impose the boundary conditions on the gravitational potentials at 
the surface of the star rather than in the asymptotic wave zone of the gravitational field. This 
approximation is tested here by predicting the frequencies of the outgoing non-radial hydrodynamic 
modes of non-rotating stars. The real parts of the frequencies are determined with an accuracy that 
is better than our knowledge of the exact frequencies (about 0.01%) except in the most relativistic 
models where it decreases to about 0.1%. The imaginary parts of the frequencies are determined 
with an accuracy of approximately M/R, where M is the mass and R is the radius of the star in 
question. 
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I. INTRODUCTION 

The pulsation of stars in general relativity theory 
has been the subject of scholarly investigation for many 
years. Consequently there are fairly well developed theo- 
ries for modeling these oscillations 0] ||^ |^ . In com- 
parison with their Newtonian counterparts, these stellar 
pulsations are fundamentally different, more interesting, 
and rather more difficult to evaluate. These differences 
arise because these pulsations couple to gravitational ra- 
diation. The illumination of a star by gravitational ra- 
diation, for example, causes that star to oscillate at the 
frequency of the incoming radiation. Thus, stars in gen- 
eral relativity theory can oscillate at any frequency at 
all! 

The outgoing (or more commonly but less descriptively, 
the quasi-normal) modes of relativistic stars are defined 
by an additional boundary condition: that they oscillate 
without incoming gravitational radiation. This boundary 
condition is simple enough to understand; however, in 
practice it is rather difficult to impose numerically when 
finding the modes of realistic relativistic stellar models. 
This difficulty arises because the external gravitational 
perturbations must be decomposed into incoming and 
outgoing parts so that solutions having only outgoing ra- 
diation can be identified. This can only be done simply in 
the region far away from the star in the wave zone of the 
gravitational field. In non-rotating stars it is straightfor- 
ward to integrate the exterior gravitational perturbations 
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out into the wave zone and perform this decomposition 
In that case the exterior gravitational field is the 
simple analytic Schwarzschild geometry, and the gravita- 
tional perturbations obey reasonably simple ordinary dif- 
ferential equations. In the case of rotating stars, however, 
this straightforward approach is not possible. The exte- 
rior geometries of rotating stars in general relativity are 
not known analytically. They are determined on numer- 
ical grids that become asymptotically poor as they ap- 
proach the wave zone. The system of partial differential 
equations that determines the exterior gravitational per- 
turbations of an oscillating star have never been solved 
numerically on such a grid. Thus, the study of the pul- 
sations of rotating stars to date has been limited to the 
Newtonian Q and post-Newtonian Q approximations. 

The purpose of this paper is to explore the possibility 
of replacing the outgoing-radiation boundary condition 
on the gravitational potentials with a condition that is 
easier to implement numerically. We derive conditions 
that are satisfied approximately by the outgoing gravi- 
tational perturbations in the near zone of the gravita- 
tional field. These conditions can easily be imposed on 
the gravitational fields in the well determined spacetime 
region just outside the star, even in numerically deter- 
mined models. We propose that these near-zone bound- 
ary conditions be imposed on the gravitational potentials 
at the surface of the star as a substitute for the outgoing- 
radiation boundary condition. We test the accuracy of 
this proposal by evaluating the non-radial modes of non- 
rotating stars, which can also be determined exactly. We 
find that the hydrodynamic modes are determined using 
the near-zone boundary condition with considerable pre- 
cision: the real parts of the frequencies are determined 
with an accuracy better than our knowledge of the exact 



1 



frequencies (about 0.01%) except in the most relativistic 
models where it decreases to about 0.1%. The imaginary 
parts are determined with an accuracy of about M/R for 
typical neutron-star equations of state. 

As part of this analysis we have derived a new repre- 
sentation of the exact relativistic pulsation equations for 
non-rotating stars, which is of considerable interest in its 
own right. We reduce the equations to a pair of second- 
order ordinary differential equations for two scalar poten- 
tials: one fluid perturbation potential, and one gravita- 
tional potential. This representation of the equations is 
exactly analogous to the Newtonian pulsation equations 

. The analysis of the near-zone boundary conditions 
described here is more straightforward using this new 
representation of the pulsation equations. 

Section |l^ describes the new representation of the rela- 
tivistic pulsation equations in terms of two scalar poten- 
tials. Section |l| presents the standard boundary condi- 
tions used to determine the outgoing modes of relativistic 
stars in terms of the two-potential formalism introduced 
here. Section IV derives approximate expressions for the 



gravitational potentials of the outgoing modes that are 
valid in the near zone of the star. These expressions 
are used as the near-zone boundary conditions. The ac- 
curacy of these boundary conditions is tested by deter- 
mining the modes of a family of typical neutron-star like 
relativistic stellar models. Section ^ concludes with a 
discussion of the feasibility of extending these near-zone 
techniques to rotating stars in general relativity theory. 
This paper also has two technical appendices which give 
in detail the new two-potential forms of the relativistic 
pulsation equations (Appendix A), and series solutions of 
these equations that are useful for imposing the bound- 
ary conditions within the star (Appendix B). 



II. NON-RADIAL PULSATIONS 

The pulsations of a relativistic stellar model are de- 
scribed by the perturbations in the spacetime metric 
5gabi and the associated perturbations in the stress- 
energy of the stellar matter ST"'^. The dynamical proper- 
ties of these quantities are determined by the perturbed 
Einstein equation. 
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and the appropriate outgoing gravitational radiation 
boundary condition at infinity. Here we consider in detail 
the non-radial pulsations of a spherical stellar model. It 
is convenient to describe the perturbations of the metric 
in such a spacetime using the Regge- Wheeler gauge, in 
which the metric and its perturbations have the following 
form, 

{gab + Sgab)dx''dx' = -e'^(l + Ho Yi„, e"^') df 
+2iHiYirne''^*dtdr + e^{l~HoYime''^*)dr^ 
+r^{l -KYi^ e"^*) {d9^ + sin^ 9 d(p^). (2) 



The perturbation functions Hq, Hi, and K, and the func- 
tions that describe the spacetime of the background equi- 
librium star, v and A, depend only on the radial coordi- 
nate r. The Yim are the standard spherical harmonics, 
and the constant lo is the frequency of the mode. The 
explicit expressions for the components of the perturbed 
Einstein tensor SG""^ that appear in Eq. (|l|) have been 
given numerous times elsewhere in terms of these fields 
0, and will not be repeated here. We use standard geo- 
metrical units in which G = c = 1. 

The perturbations in the stress-energy tensor ^T"*" are 
assumed here to be those of a perfect fluid: 

ST'''' = {Sp + 5p)u°-u'' + {p + p){du''u'' + u'^Su'') 

+Spg''''-pg^'g''%,a, (3) 

where p, p, and denote the density, pressure, and 
four velocity {u°'Ua = —1) of the unperturbed equilbrium 
model; and quantities preceeded by 5 denote their Eule- 
rian perturbations. For simplicity here we restrict our 
attention to simple barotropic perturbations, i.e., those 
where the density perturbation is proportional to the 
pressure perturbation: 



Sp=±Sp=^Sp{r)Y,,^e^-'. 
dp dp 



(4) 



The perturbed conservation law (5(VqT°'') = can be 
solved analytically to obtain a simple and useful expres- 
sion for the four- velocity perturbations. This solution is 
most conveniently expressed in terms of a scalar pertur- 
bation function SU, defined as 



SU{r) = M!i + iHo(r). 

p + p 



(5) 



Solving the spatial part of the perturbed conservation 
law gives an expression for du°' in terms of SU and the 
metric perturbations S: 



Su''= ^e^/^e^^'ViSUYir.^) 

+ (ie'^/2^oV"i + e-''^^HiV''r)Yi 



(6) 



This expression allows us to describe a stellar oscillation 
completely in terms of the four scalar functions 6U, Hq, 
Hi and K. 

The perturbed Einstein equation (|l|) can be used to re- 
duce further the number of independent functions needed 
to determine the perturbed state of the star. The 6G^g 
and SG^e components can be used to express Hi and H[ 
in terms of the other fields: 



'Hi =H'-K' + v'Hq, 



(7) 



ujH'i = -\uj{v' - X')Hi - uj^e^{Ho + K) 

-lQTT{p + p)e''+^6U. (8) 

We use the notation ' to denote differentiation with re- 
spect to r. The 6G^r and the 5G^t components together 
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with Eqs. (Q) and can similarly be used to express K 
and K' in terms of Hq, Hq, 5U and SU': 

K = aiH'a + aa-ffo + a^SU' + a^6U, (9) 



K' ^ PiH'^ + + MU' + PiSU, (10) 

where the functions and Pi depend on the equilibrium 
structure of the star (and the frequency) and are given 
explicitly in Appendix These expressions, Eqs. (|^)- 
( p^ ) , can be used to eliminate Hi, H[, K and K' from all 
of the remaining pulsation equations. The equations that 
determine the modes of relativistic stellar models can be 
reduced therefore to a pair of second-order equations for 
the fields Hq and 5U . A second-order equation for is 
obtained by using the derivative of Eq. (Q) to eliminate 
K" from the 5G*t component of Einstein's equation. A 
second-order equation for 5U is obtained from the time 
component of the perturbed stress-energy conservation 
law. In each of these equations all occurances of Hi and 
H'l are eliminated with Eqs. (R) and (0), while K and 
K' are eliminated using Eqs. (g) and (p^. The resulting 
second-order equations can be expressed in the form: 



5U"- 



+ vi SU' 



2 v' dp 
r 2 dp 
UJ^dp _ l{l + l) 



V2 



e^6U 







I" dp 

2e^dp 



- V4 



(11) 



asymptotic behavior of the perturbed fields from this rep- 
resentation of the equations. The functions Vi and 77^ are 
smaller than the terms explicitly given in Eqs. (^l]) and 
( |l^ ) by a factor of order M/R, with M the total mass 
and R the total radius of the star. Thus these terms 
are small except in the interiors of extremely relativis- 
tic models. A simple and elegant form of the relativistic 
Cowling approximation can be obtained by setting 

i/o = in Eq. (pT|). In this case the single fluid potential 
6U is determined by a single second-order wave equation. 
The solutions to this equation give fairly good approxi- 
mations of the hydrodynamic modes of relativistic stars 
pl| . Similarly we presume that approximate expressions 
for the even-parity w-modes (or 'wave' modes [Q) could 
be obtained by setting 6U = in Eq. ( |l2|) and solving 
the resulting homogeneous second-order equation for Hg . 



III. BOUNDARY CONDITIONS 



The physical solutions to Eqs. ( |ll| ) and (12) are identi- 
fied by imposing suitable boundary conditions at the cen- 
ter and surface of the fluid interior of the star, and also in 
the asymptotic wave zone of the gravitational field. The 
boundary conditions associated with the stellar fluid re- 
gion are rather simple. At the center of the star the 
perturbations are simply required to be finite. This con- 
dition eliminates one singular solution to the equations 
for each of the functions Hq and SU. The non-singular 
solutions can be approximated as power series near r = 0, 
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H 



W2 l{l+l) dp 

dp 



msu' 



e^SU 



(12) 



where the functions Vi and rji depend on the structure of 
the equilibrium star (and the frequency) and are given in 
Appendix 

Equations ( [ll| ) and (|l2|) are the general equations for 
the linear (even-parity) modes of general relativistic stel- 
lar models. This form of the equations has several ad- 
vantages over earlier representations of the fourth-order 
system needed to describe these oscillations Q , , [0 . 
Equations (|l^) and ( |l2|) reduce in a straightforward way, 
for example, to the Newtonian equations for the non- 
radial modes in the appropriate weak-field slow-motion 
limit. (Simply set Vi — rji ^ 0, = —2p'/p, e'^ = ^ 1 
in Eqs. and |l|], and ^ p = in Eq. ||l||.) The 
equation for the gravitational perturbation Hq reduces 
to a scalar wave equation in the region of spacetime far 
away from the star, and in the static limit to Laplace's 
equation. Thus it is straightforward to determine the 



Hq = Br' + 0{r 



l+2\ 



(14) 



where A and B are arbitrary constants. 

At the surface of the fluid region a boundary condi- 
tion must also be imposed that fixes the location of the 
perturbed surface of the star. Those thermodynamic po- 
tentials that vanish continuously on the surface of the 
equilibrium star (e.g. the pressure) must also vanish on 
the moving surface of the perturbed star. This condition 
is equivalent to the requirement that the Lagrangian per- 
turbation of these potentials must vanish at the surface of 
the star jl^]. The thermodynamic potential h{p) defined 
by 



Hp) 



dp' 



p{p') +p' 



(15) 



is the ideal choice to implement this boundary condition 
]l3[ . The needed condition is traditionally written Ah — 
on the surface of the star. This condition can also be 
re-expressed in terms of the Eulerian quantities as 



Ah = Sh+ - — Su^Vah = 0, 
iu) 



(16) 
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for the perturbations of a non-rotating star considered 
here. Since Sh = Sp/{p + p), this expression can be re- 
written in terms of the fields SU and Hq: 



SU' - (1 - (3M - - P2)H„ 



= 0. 



(17) 



In Appendix |^ power series expressions (vaUd near the 
surface of the star) are presented for the fields Hq and 
6U which satisfy this condition. These series solutions are 
useful for imposing the boundary condition numerically 
on a finite grid having points near, but not on the actual 
surface of the star. 

The boundary conditions associated with the fluid in- 
terior of the star, Eqs. (|l3|), jl^), and (|T7|), are sufficient 
to fix a unique (up to overall scale) solution to Eqs. ( pT| ) 
and (|l2|) for each value of the frequency ui. These con- 
ditions do not however fix the frequency of the pulsation 
modes of these stars. Indeed, relativistic stars can oscil- 
late at any frequency at all if they are driven by incoming 
gravitational radiation of that frequency. The outgoing 
modes of relativistic stars are defined by the additional 
condition that they contain no incoming gravitational ra- 
diation. This condition must be imposed as an additional 
boundary condition on the gravitational potentials in the 
asymptotic wave zone. 

In the exterior of the star the gravitational potential 
Hq is determined by the vacuum {p = p = ij^ = rj^ — 0) 
limit of Eq. (^2|). The behavior of the solutions in the 
asymptotic wave zone are of particular interest to us. 
Thus, we consider the limiting form of this equation when 
M/r <^ 1, where M is the total mass of the star. In this 
limit Eq. reduces to the homogeneous equation. 



-+Vi]Ho 



Hn=0, 



(18) 



where 771 and 772 are given in this case by 



t|(/ - 1)(/ + 2)A[l{l + 1) - 2r^uj^] - 1}, (19) 



where ji and yi are spherical Bessel functions and C and 
D are arbitrary constants. In the asymptotic wave zone, 
ro; ^ 1, the Bessel functions have simple asymptotic 
forms: 



[Cjiircj) + Dyiircj)] e*' 



D + iC 



^iuj{t-\-r) — il7T /2 



2ruj 



D-iC 
2roj ' 



^iuj{t — r)-\-il7r/2 



(23) 



The condition that the solution contain no incoming ra- 
diation reduces to the constraint D = —iC. Thus the 
required boundary condition on the gravitational poten- 
tial Hq is that it approach (as r 00) the asymptotic 
expression 

/(/-l)(2/-l)n [^ + ^(-n- (24) 
where C is an arbitrary constant. 




M/R 

FIG. 1. Frequencies of the fundamental hydrodynamic 
modes as a function of M / R for typical relativisitic stellar 
models. 



r]2 = ler^cJ^'A, 



and where A reduces to 
1 



(20) 



[(/ - 1)(; + 2) - 2r^uj^] [l{l + 1) - 2r^LU^] + Ar'^u^. 

(21) 

Remarkably, the general solution to this equation has 
been found analytically ||l4|: 



Hq 



r^-^ + l + i/(/ + 1)^(^2 

2iuj'+^[Cji{nj)+DyiirLj)] 
^ l{l-l)(2l-iy.\ 



(22) 



IV. NEAR-ZONE BOUNDARY CONDITIONS 

The frequencies of the hydrodynamic modes of inter- 
est to us here are approximately related to the mass 
M and radius R of the star by R'^uj'^ w IM/R. Fig- 
ure illustrates the frequencies of these modes for "typ- 
ical" relativistic stellar models computed from the sim- 
plified neutron-star like equation of state p = Kp^, with 
K = 1.3346 X lO^cm^/gs^. (The maximum mass model 
for this equation of state has M = I.QIIMq, and the ra- 
dius of the I.4M0 model is i? = 14.45km.) The surfaces 
of these stars are approximately in the near zone of the 
gravitational potential, r^u;^ < 1, and also in the zone 
(an extremely good approximation for the more Newto- 
nian models) where M/r <C 1. The expression for Hq 
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given in Eq. ( P4 ) is valid in the entire exterior region 
of the star where M/r <C 1. Thus, it is also valid in a 
portion of the near zone for these modes. In this region 
the expression for can be simplified by employing the 
asymptotic expansions for the spherical Bessel functions 
< 1: 



C 



-{1 - iiV,(rw)2'+i + 0{rujf + 0(A//r)}, (25) 



where Ni is given by 



N, 



;(/-l)(2; + l)[(2Z-l)!!] 



(26) 



The lowest-order version of the near-zone boundary 
condition sets i?o to the value given by Eq. (p5|): 



C 



1 - iNi{ruj) 



21+1 



(27) 



The value of the derivative Hq is obtained by differentia- 
tion of Eq. (|2^). The next question is, "Where should 
these conditions be imposed?" Moving the boundary 
away from the star improves the accuracy of the condition 
M/r <^ 1 used to derive Eq. (^), but reduces the accu- 
racy of the condition r^w^ ^ 1. As Fig. |l| illustrates, the 
r^w^ < 1 condition is just barely satisfied on the surfaces 
of the most relativistic stellar models. Thus, the most 
appropriate place to impose the near-zone boundary con- 
dition, Eq. (^, is directly on the surface of the star. We 
have solved the relativistic pulsation equations (pT|)~(|T^) 
numerically using these lowest-order near-zone boundary 
conditions for the stellar models described above. Fig- 
ures ^ and H illustrate the real and imaginary parts of 
the frequencies determined in this way. The frequencies 
in these figures are given in units of \J MjB? . Also de- 
picted in these figures are the frequencies obtained by 
solving Eqs. ([ll])-(|l^) with the exact outgoing-radiation 
boundary condition, Eq. (^^, imposed far away from the 
star in the wave zone. These figures illustrate that the 
near-zone boundary condition does reproduce in a quali- 
tative way the frequencies of the hydrodynamic modes of 
these stars. For the more Newtonian models, MjR <^ 1, 
the frequencies are determined quite accurately (as ex- 
pected). For the more relativistic models, however, the 
agreement is not as good: the real parts of the frequen- 
cies agree with the exact to within about 4%, while the 
imaginary parts only agree to within about an order of 
magnitude. This level of accuracy is, nevertheless, some- 
what better than that obtained with the standard post- 
Newtonian approximation |q]. 
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FIG. 2. Real parts of the frequencies computed with 
the lowest-order near-zone boundary condition as given in 
Eq. 
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FIG. 3. Imaginary parts of the frequencies computed with 
the lowest-order near-zone boundary condition as given in 
Eq. 



The expression for Hq given in Eq. (g5|) includes only 
the lowest-order terms in (roj)^ and M/r. By keeping 
additional terms we could in principle construct a se- 
quence of higher-order (and one might expect more ac- 
curate) near-zone boundary conditions. We have found 
(after considerable numerical experimentation) that the 
real part of Hq is most sensitive to the M/r terms in 
this expansion, while the imaginary part is most sensi- 
tive to the (rw)^ terms. This makes a certain amount 
of physical sense: the real part of Hq most strongly in- 
fluences the real part of us. This hydrodynamic aspect 
of the mode is most sensitive to the static "Newtonian" 
M/r behavior of the potential. The imaginary part of 
Hq determines, in effect, the imaginary part of uj. This 
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radiation- reaction aspect of the mode is most sensitive to 
the dynamical (rw)^ behavior of the potential. Thus, we 
have constructed a boundary condition that keeps the 
M/r terms to all orders for the real part of iJo, while 
keeping the {rtoY terms to all orders for the imaginary 
part of Hq: 



Ho = C{Q\{x) 



1 



(/-8)M 



2l{2l - 1) 



-iK, 



(28) 



where C is an arbitrary constant, x 
is given by 



K, 



2{l + 2y.{Muj) 



i+i 



i{i~i){2i~iy.\{2i + iy. 



r/M - 1, and Ki 



(29) 



The function Qf{x) is the associated Legendre function 
of the second kind. It is the exact solution to the static 
(w = 0) limit of the Hg equation, Eq. (|T^), that falls 
off like l/r'+^ for large values of r |l5|. Thus, the real 
part of Hq is given in Eq. (^8|) by an expression that is 
exact to all orders in M/r when uj — 0. The lowest-order 
(rw)^ correction has also been added to this term. The 
imaginary part of Hq in Eq. (|28|) is taken to be the exact 
solution to the Hq equation (|12[) when M = 0. This is 
just the spherical Bessel function solution that is given in 
Eq. ( p^ ) with the outgoing-radiation boundary condition. 
This solution is exact to all orders in (rw)^ when M = 0. 
The constant Ki was chosen so that the expression in 
Eq. ( p^ ) approaches (|2^) in the appropriate M/r <^ 1 
and (ruj)'^ ^ 1 limit. 

We have solved the relativistic pulsation equations 
(|ri|)^(P^ numerically using the higher-order near-zone 
boundary condition given in Eq. (Kq) imposed at the 
surface of the star r = R. The derivative of Hq is set 
to the derivative of Eq. ( [2^ ) at r = i?. Figures I and I 
illustrate the real and imaginary parts of the frequencies 
computed in this way. The real parts of the frequen- 
cies computed with the higher-order near-zone bound- 
ary conditions agree with those computed with the exact 
outgoing-radiation boundary condition with impressive 
accuracy. The fractional error is less than about 0.01%, 
the numerical accuracy with which we compute the exact 
frequencies, except in the most relativistic models. For 
the I = 2 /-modes this error is about 0.1% in the maxi- 
mum mass model, and about 0.02% for the 1.4Mq model. 
The accuracy of the real parts of the frequencies increases 
as / increases. The imaginary parts of the frequencies are 
determined with accuracy of about M/R; considerably 
better than was obtained with the lowest-order near-zone 
boundary condition. The maximum error is about 33% 
for the I — 3 mode of the most relativistic model. The 
accuracy for the modes of the IAMq models is about 
10%. This level of accuracy in the imaginary part of the 



frequency is consistent with the fact that we have ignored 
terms of order M /r in Eq. (H) . 




M/R 

FIG. 4. Real parts of the frequencies computed with 
the higher-order near-zone boundary condition as given in 
Eq. (pi). 
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FIG. 5. Imaginary parts of the frequencies computed with 
the liigher-order near-zone boundary condition as given in 
Eq. (|2|). 

To further illustrate the accuracy of the near-zone 
boundary conditions we present in Table ^ the frequen- 
cies for the modes of the 1.4Af0 models computed with 
several different approximation methods. The Cowling 
approximation used here is obtained by setting Hq = in 
Eq. (|ll|) . The post-Newtonian frequencies presented here 
are obtained from the 7 = 2 post-Newtonian frequencies 
tabulated by Cutler and Lindblom |^ and adjusted to the 
particular model considered here. The near-zone bound- 
ary conditions were used in both the lowest-order form 
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of Eq. (H^) and the higher-order form of Eq. (|2|) . 

The detailed numerical results presented here were ob- 
tained for the 7 = 2 polytropic equation of state de- 
scribed above. We have also evaluated the modes of neu- 
tron star models based on a softer 7 = 5/3 polytropic 
equation of state, and also on the more realistic Bethe- 
Johnson equation of state ||l8| . We find that the accuracy 
of the frequencies in these two additional cases is com- 
parable to (and in fact slightly better than) the 7 = 2 
polytropic case discussed in detail here. 

V. CONCLUDING REMARKS 

This paper presents a new formulation of the relativis- 
tic pulsation equations in terms of two scalar potentials: 
one fluid and one gravitational perturbation field. These 
potentials satisfy a pair of second-order ordinary differen- 
tial equations that are analogous to the Newtonian pulsa- 
tion equations. These equations are quite general: they 
describe the complete set of even-parity modes of rela- 
tivistic stars, including the w-modes. This formulation 
of the relativistic pulsation equations has a significantly 
simpler structure than earlier representations. 

The primary purpose of this paper was to explore the 
possibility of replacing the outgoing-radiation boundary 
condition imposed in the wave zone with appropriate con- 
ditions imposed on the surface of the star. We have used 
the two-potential forms for the stellar pulsation equa- 
tions to derive the appropriate near-zone boundary con- 
ditions on the gravitational potential. We have shown 
that the outgoing hydrodynamic modes can be deter- 
mined with considerable precision using these near-zone 
boundary conditions. Unfortunately this success is some- 
what vacuous: we already knew how to evaluate these 
modes exactly. However we believe that this proceedure 
should also allow us to evaluate the modes of rapidly ro- 
tating fully relativistic stars with a precision similar to 
that obtained here. The near-zone boundary conditions 
that we present here depend only weakly on the struc- 
ture of the exterior gravitational field of the star. The 
boundary condition used to determine the real part of 
the potential Hq does depend on the exact static solu- 
tion to the exterior gravitational perturbation equations. 
But such static solutions have been found numerically in 
rapidly- rotating strongly relativistic models . We be- 
lieve that these numerical solutions can be used to formu- 
late a near-zone boundary condition that will be accurate 
even for rapidly rotating models. The boundary condi- 
tion on the imaginary part of Hq does not depend on the 
structure of the external gravitational field of the star at 
all: it is just the flat-space perturbation representing an 
outgoing wave. Finally, the modes of primary interest in 
rotating stars have the property that their frequencies lu 
decrease toward zero as the angular velocity of the star 
increases. For these modes the near-zone approximation 
r^o;^ ^ 1 on which our boundary condition is based will 



be satisfied more exactly as the angular velocity of the 
star increases. 
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APPENDIX A: THE PULSATION EQUATIONS 



In Section |l| are presented the general forms of the equations that express K and K' in terms of iJo, H'^, SU and 
SU': 



K' = /3ii7^ + (32Ho + MU' + PiSU. 
Here we present the complete expressions for the coefficients and f3i: 



ai = -rA{4r2w2e-'^ + rv' [IQirr'^ {p + p) - l{l + 1)]}, 
a2 = A{r(2 - rv')[l&T:r'^ [p + p)i^' - l{l + l)t/' + iru"^ 

~[l{l + l)e^ - 2(1 - rv') - Snr^ip + p)e^][16nr^{p + p) - 1(1 + 1) + 2r^uj'^e-'']} , 
as = -167rr^A(p +p)(2 - riy'), 

a4 = -16TTr^{p + p)e^A[l6Trr^{p + p) - l{l + 1) + 2r^uj^e-''], 
/3i 

Pi 



(Al) 
(A2) 

(A3) 

(A4) 
(A5) 
(A6) 
(A7) 



A{[2r^Lj^e-'' - {I - + 2)][167Tr\p + p) - l{l + l)]e^ + 2r^Lj^e-'' {2 ~ ri^')} , 
A{rLu^e-''{2 - rv')[l{l + l)e^ - 2(1 - ru') - 8nr^{p + p)e^] 

-[{I -l)(l + 2)- 2r^uj^e"'][l&nr^v'{p + p) - l{l + + 2rLu^e-'']e^) , 
167rr2(p + p)Ae^[(; - 1)(; + 2) - 2r^uj^e-% 
16nr^uj'^{2 - riy'){p + p)Ae^-'' , 

where A is defined by 

A-i = [{I - 1)(/ + 2) - 2r^iu^e-''][l{l + 1) - 2r^iu^e-'' ~ IGnr^ip + p)]e^ + r^Lu^e-''{2 - ri^'f 

This transformation becomes singular whenever A~^ vanishes. We do not yet understand the general conditions 
under which this singularity can occur (if any). We have determined numerically, however, that A > for all r in the 
hydrodynamic modes studied here. 

Also in Section || were presented the general forms of the second-order equations for 6U and Hq: 



(A8) 
(A9) 
(AlO) 

(All) 



5U" 



v' dp 



2 dp 



+ vi SU'+ 



V2 



3^SU 



bj^ dp 
2e'' dp 



Hi 



LU^ 1(1 + 1) ^ . .dp 
— - , +47r(p+p)-f +7^2 



= m5u' - 

Here we give the complete expressions for the coefficients Vi and rji 

"2 v'fdp 



dp 

dp 

8n{p + p} — +r]i 
dp 



e^6U. 



V2 = 16tt{p + p) + e 



r 2 

A 



1^3 



I/' / dp 



2 \dp 



dp 

2 v'fdp 
r '2\dp 

{I -Pi), 

v' ( dp 

dp 



P3, 



2 

2 



1 



Pi, 



- P2), 



-\') + }i2-ru'){P,-l), 
= e-^{v" + (v')^ + ^(2 - rv')[4P2 - v' + X']} + An{p + p), 
m = -l{2-riy')P3, 

m 



-8TT{p + p)-U2~rv')e-^P4 



(A12) 

(A13) 

(AM) 
(A15) 
(A16) 

(A17) 

(A18) 
(A19) 
(A20) 
(A21) 
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APPENDIX B: POWER SERIES SOLUTIONS 



SU{r) = SU{R) - (i? - r)SU'{R) + 0{R - r)i+^, (B8) 



It is helpful to have power series solutions to Eqs. (|l l| ) 
and (|l^) for the functions 5U and that satisfy the 
appropriate boundary conditions near the center, r = 0, 
and near the surface, r = i?, of the star. Such series 
solutions are useful, for example, in allowing the bound- 
ary conditions to be imposed numerically at grid points 
near, but not on the actual boundaries. Near r = the 
situation is straightforward: the functions 5U and Hq are 
given by the series solutions 



5U = Ar' +0(r'+2). 



Hn = Br' 



0{: 



J+2\ 



(Bl) 



(B2) 



where A and B are arbitrary constants. 

Near the surface of the star, however, the situation 
is more delicate. The thermodynamic derivative dp/ dp 
typically diverges at the surface of the star. Thus the 
second derivatives of SU and Hq may be infinite at r = 
R. This makes a straightforward power series expansion 
impossible. For definiteness we assume here that in the 
neighborhood of the surface of the star the equation of 
state can be represented as a polytrope 



(B3) 



where k and 7 are constants. Near the surface of such a 
star the equilibrium structure equations can be solved as 
power series expansions. From this we learn that 



dp _ p _ R{R - 2M) 
dp~ rP~ (7 - l)M{R - r) 



[l + 0{R-r)], (B4) 



where < t < 1 is a positive constant that depends on 
the equation of state. For soft equations of state 1 < 7 < 
2 the singularities are relatively beneign, and the second 
derivatives H'^ and 5U" are finite at r ~ i?. Thus, in this 
case series expansions analogous to Eqs. (B7) and (^^ ) 
can be obtained for the derivatives: 



H'^{r) = H'^{R) - (i? - r)H'^{R) + 0{R 



(B9) 



5U'{r) = 5U'{R) ~{R- r)5U"{R) + 0{R - rf, (BIO) 

where the constants HJUR) and SU"{R) are evaluated 
directly from the Eqs. ( pT^ and ([T^). In the stiff case, 7 > 
2, series expansions can also be found using Eqs. (B4)- 
@ and the equations These expansions are 

not analytic at r = i? but can nevertheless be used to 
approximate the derivatives HQ(r) and 5U'(r) near the 
surface: 



H'„{r)= H',{R) 



^R^p{r)[2SU{R)^Ho{R)] 
x[l + 0{R-r)], 



(Bll) 



SU'{r) = SU'{R) 

-i{|f _ f3,(^R)]p(r)[2SUiR) - Horn 

+l3z{r)5U'{R) + l3^{r)6U{R)Yl + 0{R - r)]. (B12) 



7 - 1 M{R - r) 



7K R{R-2M) 



[l + 0{R~r)]. (B5) 



Thus the derivative dp /dp that appears as a coefficient in 
Eq. always diverges for polytropcs. The combination 
{p + p)dp/dp that appears in Eq. ( |l2[ ) diverges for stiff 
equations of state (7 > 2) but not for softer equations of 
state: 



iP + p) 



1 



dp 
dp 



7 - 1 M{R-r) 
7K R{R-2M) 



[l + 0{R-r)]. 

(B6) 



The boundary conditions determine the surface values 
of 5U'{R) and H'^{R). The value of H'^{R) is determined 
by continuity at r = i? from the outgoing solution in the 
exterior of the star. The value of dU'{R) is determined by 
the boundary condition Eq. ( p^ . These constants can be 
used therefore to give linear expressions for the functions 
Ho{r) and 6U{r) near the surface: 

Ho{r) = Ho{R) -{R- r)H'^{R) + 0{R - rf+\ (B7) 
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TABLE I. The exact frequencies of the /-modes of IAMq stellar models are compared with 
those computed using the near-zone boundary condition, and with those computed using several 
other approximation methods. 





I 


Cowling 


post-Newtonian 


Lowest-Order 
Near-Zone 


Higher-Order 
Near-Zone 


Exact 


Re{u;) 


2 


1.4484 


1.231 


1.1672 


1.2018 


1.2016 


Re(tj) 


3 


1.7469 


1.619 


1.5716 


1.5863 


1.5862 


Re(a;) 


4 


1.9906 


1.907 


1.8652 


1.8739 


1.8738 


Ke{uj) 


5 


2.2031 


2.147 


2.1071 


2.1130 


2.1129 


Im(a;) 


2 






5.60 X 10"* 


3.53 X lO""* 


3.95 X lO"'' 


Im(ct;) 


3 






7.29 X lO"'^ 


3.56 X 10"*^ 


3.96 X lO"'' 


Im(ct;) 


4 






9.36 X 10"** 


3.78 X 10"* 


4.01 X 10"** 


Im(ct;) 


5 






1.09 X 10"^ 


3.70 X 10"^° 


3.63 X 10"" 
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